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ABSTRACT 

We present the results of an XMM- Newton mosaic covering the central ~ 200 kpc 
of the nearby Virgo cluster. We focus on a strong surface brightness discontinuity 
in the outskirts of the brightest cluster galaxy, M87. Using both XMM-Newton and 
Suzaku, we derive accurate temperature and metallicity profiles across this feature 
and show that it is a cold front probably due to sloshing of the Virgo ICM. It is also 
associated with a discontinuity in the chemical composition. The gas in the inner, 
bright region of the front is ~40% more abundant in Fe than the gas outside the front, 
suggesting the important role of sloshing in transporting metals through the ICM. For 
the first time, we provide a quantitative estimate of the mass of Fe transported by a 
cold front. This amounts to ^6% of the total Fe mass within the radial range affected 
by sloshing, significantly more than the amount of metals transported by the AGN 
in the same cluster core. The very low Fe abundance of only ~0.2 solar immediately 
outside the cold front at a radius of 90 kpc suggests we are witnessing first-hand the 
transport of higher metallicity gas into a pristine region, whose abundance is typical 
of the cluster outskirts. The Mg/Fe and O/Fe abundance ratios remain approximately 
constant over the entire radial range between the centre of M87 and the faint side of 
the cold front, which requires the presence of a centrally peaked distribution not only 
for Fe but also for core-collapse type supernova products. This peak may stem from 
the star formation triggered as the BCG assembled during the protocluster phase. 

Key words: X-rays: galaxies: clusters; galaxies: individual: M87; cooling flows; ISM: 
abundances 



1 INTRODUCTION 

In the last decades, our understanding of the intricate phe- 
nomena in clusters of galaxies has advanced tremendously 
based on observations with the constantly improving gener- 
ations of X-ray satellites. Clusters are now known to show 
complex morphologies and exciting physics on all scales and 
are thus among the most rewarding objects to target in 
the X-ray domain. Detailed spatial and spectral information 
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about clusters can reveal important clues about the ener- 
getics of large-scale gas motions in the intracluster medium 
(ICM) and about the interaction between the ICM and the 
active galactic nucleus (AGN) in the central galaxy. Both 
AGN-ICM interactions in the form of shocks and bubbles 



(Binney & Tabor 


1995 


Churazov et al. 


2000; Bohringer 


et al.|2002||Bruggen & Kaiser .2002 Birzan et al.|2004||For- 


man et al."1|2005| |Voit & Donahue||2005| |Fabian et al.||2006| 


McNamara & Nulsen||2007 


1 and gas "sloshing" in the clus- 


ter dark matter potential ( Markevitch et al. 


2001 Churazov 


et al. 2003 Markevitch & Vikhlinin 2007 


1 have been in- 



voked, among other explanations, as possible ways to heat 
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the gas in the centre of cool-core clusters, preventing the 
radiative cooling of the gas at high rates which should occur 



staggered placement of these two features is a typical exam- 



in the absence of heating but which is not observed (Pe- 



terson et al. 2001 2003 Bohringer et al. 20011. Moreover, 



detailed spatial and spectral information about clusters is 
essential in constraining the chemical enrichment history of 
the ICM (|de Plaa et al.|2007| [Werner et al.|2008||Bohringer 



fc Werner||2009[ ) and possible mechanisms responsible for 
the transport of heavy elements throughout the hot cluster 
medium ( |Rebusco et al.|2006| |Simionescu et al.|2008] |2009 1 . 

One of the best targets to perform in-depth studies of 
the ICM is M87 at the centre of the Virgo cluster. This is 
the nearest galaxy cluster (~ 16 Mpc), allowing us to resolve 
phenomena on small scales. M87 is among the brightest ex- 
tragalactic X-ray sources in the sky, a guarantee for excel- 
lent spectral statistics and accurate temperature and metal 
abundance determinations within the technical capabilities 
of the detectors. M87 is among the best studied galaxies 
in the Universe and its central region has been the tar- 
get of very deep observations with XMM-Newton, Chandra, 
and recently Suzaku. The hot gas atmosphere of M87 shows 
striking signs of AGN-ICM interaction, including an AGN- 
driven classical shock ( |Forman et al . 2005 , 2007 . Simionescu| 
et al. |2007" I and clear enhancements in X-ray surface bright- 
ness associated with the radio lobes to the east and south- 



west of the core ( jFeigelson et al.|l9 87 Bohri nger et al.|1995 
Belsole et al.||2001| | Young et al ||2002| |Forman et al. ||200~ 



2007). Initial XMM-Newton observations showed that multi- 



temperature models including a cool gas component are 
needed to describe the spectra in these regions, also known 
as the E and SW X-ray arms (|Belsole et aL|2001| |Molendi 



2002). More recently, with deeper XMM-Newton data, a cor- 
relation was found between the percentage of cool gas and 
the metallicity in these regions, based on which [Simionescu] 
|et al. | ( |2008[ ) concluded that the cool gas is metal-rich and 
that it was uplifted by the AGN (following the scenario pro- 



posed by Churazov et al. 2001 1 from the central parts of 



the galaxy which had been enriched by stellar mass loss and 
supernova explosions. The energetics and metal-transport 
mechanisms associated with the AGN activity in the central 
parts of M87 have therefore been thoroughly investigated 
to date and have yielded important information about the 
physics of the AGN-ICM interaction. Here, we present the 
results of Suzaku and XMM-Newton observations extend- 
ing the in-depth analysis to the outskirts of M87 to look for 
evidence for large-scale gas motions in the Virgo core and 
to bring gas sloshing and AGN-ICM interaction together in 
explaining the metal transport processes in the ICM. 

We focus on a clearly observed surface brightness edge 
that lies 90 kpc (19') to the north of M87 and extends 140° in 
azimuth in the archival ROSAT ~30 ks image, its large scale 
making it an optimal target for Suzaku given the limitations 
related to its point-spread function. We show, using recent 
XMM-Newton and Suzaku data, that this edge is a "cold 
front", a contact discontinuity associated with gas sloshing 
in the ICM, typical of those found in many cool core clus- 
ters ( |Markevitch fc Vikhlinin|2007[ |. [Simionescu et al.| ( |2007 l 
have previously detected a weaker counterpart cold-front 
~ 33 kpc to the south-east of the core. The opposite and 



pie of sloshing as seen in simulations (Tittley & Henriksen 
2005[|Ascasibar fc Markevitch|2006 1. As the gas sloshes back 



and forth in the gravitational potential, cooler and denser 
parcels of gas from central parts of the cluster can move and 
come into contact with the hotter outskirts creating "cold- 
front" signatures, where the temperature on the denser side 
of the surface brightness edge is lower, which is the oppo- 



site of what one expects from a shock front (Markevitch 



et al. 2001 1. Since heavy elements in cooling flow clusters are 



concentrated toward the centre (e.g. De Grandi & Molendi 
2001 Leccardi & Molendi 2008 ), the displaced parcels of cool 



central gas are also more metal-abundant. Thus, the abun- 
dance should be discontinuous across these fronts, as long as 
sloshing occurs within the region with a strong gradient. In 
general, sloshing should spread the heavy elements from the 
centre outwards, providing, beside the central AGN, a mech- 
anism for distributing and transporting the metals produced 
in the central galaxy into the ICM. 

In this paper, we discuss the abundance gradients for 
different heavy elements across the cold front to quantify the 
metal transport due to sloshing and to study the chemical 
enrichment histories of the gas on either side of the surface 
brightness edge. 



2 DATA REDUCTION 
2.1 Suzaku 

On 2008 June 8, Suzaku observed two fields offset to the 
North by 19' and 33' with respect to the centre of M87. The 
effective exposure times were 59.7 ks for the pointing fur- 
thest away from the M87 centre ( "M87Nfaint" ) and 21.9 ks 
for the pointing with a smaller radial offset ( "M87Nbright" ) . 
The distance between the Southern edge of the field-of-view 
of the M87Nfaint pointing and the surface brightness dis- 
continuity was chosen to be ~ 3' (of the order of the Suzaku 
PSF) in order to avoid contamination from the bright side 
of the front into the faint-side spectra. 

The data were reduced using the tools available in the 
HEAsoft package (version 6.7) to create a set of cleaned 
event lists with hot or flickering pixels removed and includ- 
ing only ASCA grades 0,2,3,4,6. The recommended filter- 
ing criteria were used to remove South Atlantic Anomaly 
(SAA) transits and times of low geomagnetic cut-off rigidity 
(COR), eliminating data with COR ^ 6 GV. In addition, 
all data taken within 20° of the sun-lit Earth limb or 5° 
of the dark Earth limb were excised to reduce contamina- 
tion from scattered solar X-ray flux. Data obtained with 3x3 
and 5x5 editing modes were merged into a single event file. 
The lightcurves for both data sets after the initial cleaning 
do not show anomalous count rates either in the full or in 
the soft energy bands. Ray-tracing simulations based on the 
available ROSAT data show that the expected contamina- 
tion by out-of-field-of-view emission is below 1.5% in both 
pointings. An exposure-corrected flux image of the two com- 
bined Suzaku pointings is shown in Fig. [I] (no background 
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Figure 1. Mosaic of the two Suzaku fields (exposure-corrected). 
The dashed line marks the approximate position of the outer cold 
front. 



subtraction has been performed for the image processing). 
The cold front is clearly visible in the image. 

For the spectral analysis, we subtracted the non X-ray 
background (NXB) obtained with xisnxbgen and then mod- 
eled the cosmic X-ray background (CXB) with 3 compo- 
nents: two thermal components to account for the local hot 
bubble (LHB) and for the Galactic halo (GH) emission, and 
a power-law to account for the integrated emission of un- 
resolved point sources. The temperature of the LHB was 
frozen to 0.08 keV, that of the GH to 0.2 keV, and the 
power-law index to 1.41; unless otherwise noted, the spec- 
trum normalisations were fixed based on the fluxes given by 



Kuntz & Snowden (20001 for the thermal components and 



by De Luca & Molendi 



2004) for the power-law. We note 



that the region investigated is within approximately 0.2 of 
the virial radius of the Virgo cluster, thus the CXB level is 
still low compared to the emission from the ICM. 



2.2 XMM-Newton 

A series of 6 pointings with an average exposure time of 
20 ks was performed with XMM-Newton in June 2008, with 
an additional seventh pointing of the same length completed 



on 2008 November 12. These seven pointings cover an annu- 
lus around the existing deep observations of the central part 
of M87 from 2000 June 19 and 2005 January 10, which had 
a combined effective exposure time of 120 ks. A combined 
exposure-corrected flux image of the MOS data from all the 
pointings is shown in Fig. [2] together with the residuals ob- 
tained by dividing this image by a circularly symmetric beta 
model, with a best-fit core radius of 12.6" and j3 = 0.38. 
Two surface brightness discontinuities, one at a radius of 
~33 kpc towards the SE and one at ~90 kpc towards the 
NW are clearly seen. 

For data reduction we used the 8.0.0 version of the 
XMM-Newton Science Analysis System (SAS). Out-of-time 
events were subtracted from the EPIC/pn data using the 
standard SAS prescription for the extended full frame mode. 
To remove flaring from soft protons, we used a two-stage 
cleaning of the event files. First, we extracted a light-curve 
in the hard energy band (10-12 keV for MOS1 and MOS2 
and 12-14 keV for pn) using 100 second time bins and ex- 
cluded the time periods in each observation when the count 
rate exceeded the mean by more than 3a. After this first 
cleaning, we used the resulting event files to extract a light- 
curve in a broad energy band (0.3-10 keV) with 10-s time 
bins and again excluded the time intervals with count rates 
larger than the average by more than 3a. We then examined 
the soft band (0.3-0.5 keV) images of each MOS observation 
to look for chips with anomalously high flux in this energy 
range ( jSnowden et al. 2008), and found that none of our 
observations were affected by this instrumental problem. 

For the background subtraction, we used a combination 
of blank-sky maps from which point sources have been ex- 
cised ( Read &: Ponman|2003[ ) and closed-filter observations. 
This is necessary because the instrumental background level 
of XMM-Newton is variable and increases with time. For 
each detector, we added to the corresponding blank sky 
background set a fraction of the closed filter data designed 
to compensate for the difference between the out of field of 
view (OoFoV) hard-band count rate in the observation and 



in the blank sky data (for more details, see Simionescu et al. 
[2009] ). 



3 TEMPERATURE AND IRON ABUNDANCE 
PROFILES 

We first present radial profiles across and away from the 
northwestern cold front in order to establish the projected 
temperature and Fe abundance jumps associated with this 
feature and compare the different trends towards the NW 
and SE from the core. 

For the XMM-Newton data, the spectral extraction re- 
gions used were 90°-wide sectors of elliptical annuli towards 
the NW and SE, respectively. The thickness of these annuli 
increases logarithmically with distance from the M87 centre, 
in order to balance the decreasing surface brightness. The 
exact regions used are shown in Fig. [2] For the Suzaku data, 
we chose elliptical annuli with the same ellipticity and po- 
sition angle as for the XMM-Newton data. For illustration, 
the ellipse delimiting the two Suzaku annuli immediately 



© 2009 RAS, MNRAS 000,[l}{TT] 



4 A. Simionescu et al. 




5E-05 0.0001 0.00015 0.0002 0.00025 0.00 

Figure 2. XMM-MOS mosaic of M87. Left: Exposure corrected image divided by the best-fitting radially symmetric beta model. Right: 
Extraction regions used for the radial profiles in Sect.[3]and the two Suzaku fields overplotted on the exposure corrected image. Colorbar 
units are counts per second. 



outside and immediately inside the outer cold front is over- 
plotted in Fig. [I] No sector selection was performed since the 
opening angle of the area covered by Suzaku is already lim- 
ited (see green and yellow squares overplotted in Fig. [5| and 
all annuli were fixed to have a width of 3', due to constraints 
from the Suzaku PSF. 

We use XSPEC 12.5.1 to model our spectra with a 
single-temperature plasma model in collisionally ionised 
equilibrium. Unless otherwise stated, we used the apec 
plasma model ( Smith et al.|200l l. The determination of the 
Galactic absorption column density based on X-ray spectra 
can be easily biased by subtle underlying multi-temperature 
structure in the ICM ( Simionescu et aL]|2009 1 or by uncer- 
tainties in the calibration of the oxygen edge around 0.5 keV 
and of Suzaku's contamination layer. We t herefore opt to fix 
this parameter to A h = 2.0 x 10 20 cm" 2 \hieu et al.|l996 
|Kalberla et al.|2005| >. 

Point sources identified using the X-ray images were 
excluded from the spectral analysis for the XMM-Newton 
spectra. No unusually bright point sources are present within 
the field of view of Suzaku. The spectra obtained by differ- 
ent detectors of the same satellite were fit simultaneously 
with their relative normalisations left as free parameters. 
This is mainly done in order to account for effective area 
variations due to XMM-Newton chip gaps, the failed CCD6 
on MOS1, and the different locations of the field of view 



edges for MOS and pn. The instrument normalisations of 
the three XIS detectors agree to within 5% for each region 
analyzed. The MOS, pn, and XIS1 data were fit in the 0.4- 
7.0 keV band while for XIS0 and XIS3 we used the 0.4-10 
keV band. The spectra were binned with a minimum of 30 
counts per bin before fitting. The normalisation, tempera- 
ture, and the abundances of O, Ne, Mg, Si, S, Ar, Ca, Fe, and 
Ni were left as free parameters. The abundances of the other 
elements heavier than He were fixed to 0.5 of the Solar value. 
For consistency with previous chemical enrichment studies 



of M87 with XMM-Newton by Matsushita et al. (|2003|>, we 



adopt the solar abundances of Fcldrnan (1992) 

Our spectral fitting results are shown in Fig. [3] We find 
an excellent agreement between XMM-Newton and Suzaku 
in determining all the spectral parameters, despite the differ- 
ent background subtraction techniques employed. The X-ray 
surface brightness profile, obtained by dividing the spectrum 
normalisation from XSPEC by the area of the correspond- 
ing extraction region, shows a very clear discontinuity at 
19' (~90 kpc), the radius of the outer cold front. At the 
radius of the inner cold front, the surface brightness pro- 
file towards the SE crosses from being systematically higher 
to being systematically lower than the corresponding pro- 
file towards the NW. The jump however is less pronounced, 
and is only evident once the strong radial trend is removed 
using a symmetric model, as was done in Fig. [2] The tem- 
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perature profile towards the SE rises beyond the inner cold 
front at ~33 kpc discussed by Simionescu et al. (20071 and 



then becomes flat, while the projected temperature to the 
NW stays flat or decreases slightly in the radial range be- 
yond 30 kpc and only shows a sharp rise at the cold front 
at 19' (~90 kpc). Therefore, from the temperature profiles, 
we can conclude that the surface brightness discontinuities 
both at 33 and 90 kpc are cold fronts rather than shocks, 
exhibiting a lower temperature of the gas on the denser side. 
After increasing outside of the 90 kpc cold front, the temper- 
ature seems to decrease again beyond 30'. This is however 
most likely a feature of the general Virgo temperature pro- 
file rather than being associated with the event that caused 
the cold front. A more detailed deprojection analysis will be 
presented by Forman et al. (in prep). 

Typically two types of cold fronts have been discussed in 
the literature (see Owers ct al. 2009 . Markcvitch & Vikhlinin 
2007 1 and two different mechanisms for producing them have 



been proposed. On one hand, cold fronts have been seen in 
merging clusters, where often the remnant core belonging to 
an infalling group or smaller cluster can still be identified. 
The best and most famous example here is the bullet clus- 



ter (Markevitch et al. 20021. The cold front then exists at 



the interface between the cool, dense remnant core and the 
ICM of the cluster with which it is merging. This is clearly 
not the case for M87, which appears relatively relaxed apart 
from the X-ray arms and other small scale features in the 
very center (inner ~ 1'), all of which seem related to AGN- 
ICM interaction (e.g.|Feigelson et al.|1987 Matsushita et al 



2002||Forman et al.|2 005 Sim ionescu et ai.|2008| >. The most 



probable explanation for the opposite and radially staggered 
cold fronts in M87 thus is gas sloshing, similar to that pro- 



posed by Markevitch et al. (2001 1 in A1795, where the grav- 



itational potential of a cool core cluster can be perturbed 
by a past subcluster infall. Due to this disturbance, the cen- 
tral ICM is displaced and fronts are created where cooler 
and denser parcels of gas from the centre come into con- 
tact with the hotter outskirts. Numerical simulations show 
that the sloshing typically follows a spiral pattern, which 
can explain the placement of the two cold fronts in M87 at 



different radii on opposite sides from the centre (Tittley & 



Henriksen 2005 Ascasibar & Markevitch 20061. Using the 



mass profile of Matsushita et al. (20021, the Keplerian ve- 



locity is approximately 550 km/s for an orbit with a radius 
of 33 kpc and 850 km/s for a radius of 90 kpc. A rough 
estimate for the time needed for the two fronts to become 
out of phase by 180 degrees is thus 71-/(1^33 —ujgo) = 0.4 Gyr, 
where uj — v/r is the angular velocity at the corresponding 
radius. This estimate is likely to be low. 

An interesting topic for further simulations in the case 
of M87 will be whether the large scale bulk motions due to 
ICM sloshing are strong enough to cause the observed bend- 
ing of the AGN-inflated radio lobes and of the associated E 
and SW X-ray arms. 

Apart from the marked increase in temperature beyond 
the surface brightness discontinuities at 33 and 90 kpc, an- 
other common feature is a sharp drop in Fe abundance out- 
side each cold front. A discontinuity in the metal abundance 
profile is expected because the Fe abundance in Virgo, as in 



the vast majority of cooling core clusters, is centrally peaked 
(Fig. [3] see also e.g. |Leccardi fe Molendi|2008[ ). This means 
that the central gas displaced by sloshing is also richer in 
metals and thus sloshing contributes to transporting the 
heavy elements from the centre outwards. 

To quantify this, we calculate the excess mass of Fe to- 
wards the SE behind the inner cold front (in an opening 
angle of 90° between ~14 and 33 kpc) and the excess mass 
of Fe towards the NW behind the outer cold front (in an 
opening angle of 90° between ~33 and 90 kpc) . This is done 
simply by estimating the gas mass in each annulus based 
on the spectrum normalisation, multiplying by the Fe mass 
fraction determined from the Fe abundance, subtracting the 
corresponding value for the opposite sector, and finally sum- 
ming over all annuli in the cited radial ranges. Note that this 
is likely an underestimate of the mass of Fe transported by 
the outer cold front. This is firstly because this front extends 
over an opening angle larger than 90° (although it is less 
pronounced outside the region considered here). Secondly, 
the low Fe abundance outside the outer cold front suggests 
that the initial Fe profile to the NW before the metal trans- 
port was lower than estimated from taking the SE profile at 
corresponding radii. 

We obtain an excess Fe mass of 4.5 x 10 6 Mq behind 
the inner cold front and 13.7 x 10 6 Mq behind the outer cold 
front. This represents 31% and 22% of the average Fe mass 
per sector within the corresponding radial range for the in- 
ner (33 kpc) and outer (90 kpc) fronts, respectively. If the 
Fe mass averaged between the NW and SE sectors is repre- 
sentative also for each of the other two sectors towards the 
NE and SW, then the mass of Fe transported by sloshing 
represents ~6-8% of the total Fe mass within the consid- 
ered radial ranges. Alternatively, if we relate the mass of Fe 
transported by the two cold fronts to the total Fe masses 
between 14-33 and 33-90 kpc obtained using the cumula- 
tive Fe mass profile of M87 from ( Bohringer et al.|2004 1 , we 
obtain very similar ratios of ~6-9%. 

It is noteworthy that both the large-scale 90 kpc cold 
front and the inner cold front seem to transport a very 
similar fraction of the mass of Fe available within the ra- 
dial region affected by the sloshing. By comparison, the Fe 
mass uplifted by the AGN in the bri ght X-ray arms is only 
1.5 x 10 6 M Q ( |Simionescu et al.|200"8| . Note that this is only 
an estimate of the current amount of metals being displaced 
by either process - depending on the strength of the mixing 
within the ICM, a fraction of these metals may eventually 
fall back to their original position within the cluster. 

While it was expected that the sloshing that causes cold 
fronts is an important means of transporting metals in the 
ICM and that this would cause the Fe abundance to be dis- 
continuous at the surface brightness jump towards the NW, 
the strength of this feature is surprising. Immediately out- 
side the cold front at 90 kpc, we reach a very low, flat abun- 
dance plateau. The average abundance here is 0.25 ± 0.02 
solar (sec also next section) relative to the solar units of 
Feldman (19921, which translates to only 0.17 ± 0.01 solar 



in the units of Anders & Grevesse (19891. This is as low as 



the expected metallicity close to the virial radius in other 



clusters (0.2 solar, Fujita et al. 20081. Such low Fe abun- 
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dances are rarely measured in the ICM. For a sample of 
nearby galaxy clusters observed with XMM-Newton, |Lecca-| 
rdi & Molendi ( 2008 ) find that the average metallicity profile 
keeps decreasing out to very large radii, reaching ~0.2 solar 
only at around 0.3-0.4 riso. A similar trend is shown by the 



Fe abundance profiles measured with Suzaku (e.g. Sato et al 



20091. For an average temperature of 2.5 keV, the virial ra- 



dius of the Virgo cluster based on the scaling relations of 
e.g. Arnaud et al.| ([2005) is r^oo ~ 1.2 Mpc, which makes it 
very surprising to see gas with such a low metallicity out- 
side a radius of only 90 kpc. Towards the SE, the metallic- 
ities at corresponding radii beyond the outer cold front are 
systematically higher, closer to what one would expect at 
these relatively small distances from the cluster centre. The 
likely scenario is that sloshing in the more distant past has 
influenced the abundances of the gas towards the SE and, 
because Virgo is a relatively young and unrelaxed object, 
we are only now witnessing the dispersion of higher metal- 
licity gas into the pristine region towards the NW, whose 
abundance is that typical of the cluster outskirts. 



4 METAL ABUNDANCE RATIOS AND 
SUPERNOVA ENRICHMENT HISTORY 

It is believed that the peak in Fe abundance associated with 
cool core clusters is produced by type la supernova ejecta 
in the brightest cluster galaxy (BCG) which is spatially co- 



incident with this peak (De Grandi et al. 2004 Bohringer 



et al.||2004| . While SN la have a longer delay time and still 
occur in these BCGs, which have typically old stellar popula- 
tions, core-collapse supernovae (SNcc) usually explode only 
within a short time after a star formation event and would 
thus not be expected to enhance abundances significantly 
in cool cores. Elements which are predominantly produced 
by SNcc, such as O, Ne and Mg, or even intermediate ele- 
ments produced by both SNcc and SN la such as Si and S, 
are therefore expected to be less abundant compared to Fe 
in the vicinity of the BCG. Whether or not such trends in 
the e.g. O/Fe or Si/Fe ratios are in fact observed is, how- 
ever, a matter of debate in recent literature. |Finoguenov| 
et al. ( 2000 1 show an increase of Si/Fe with radius in a sam- 
ple of clusters observed with ASCA. While the presence of 
such a gradient was confirmed, more recently, for the case 
of galaxy groups ( Rasmussen fe Ponman|2009 1, the constant 
Si/Fe and increasing O/Fe between the centre and outskirts 



of M87 prompted Finoguenov et al. ( 2002 I to propose a di- 
versity of SN la: faint SN la with less complete burning and 
higher Si/Fe yields would contribute to the enrichment of the 
central parts of the ICM, while brighter SN la with lower 
Si/Fe yields would dominate the enrichment in the outskirts. 



A strong increase in O/Fe with radius claimed by Tamura 



et al. 



( 2004 1 has not been confirmed by recent analyses of 
deeper cluster observations, which show no significant in- 
crease in O/Fe within 0.1 r2oo ( Simionescu et al.|2009 l. The 
latest Suzaku results also indicate constant Mg/Fe, Si/Fe, 
and S/Fe abundance ratios with radius and only a very weak, 
low-significance increase in O/Fe out to as far as 0.4 riso for 
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Figure 3. X-ray surface brightness, temperature and Fe abun- 
dance profiles towards the NW and SE as a function of the average 
semi-major axis (SMA) of the elliptical annuli used as extraction 
regions. The positions of the two cold fronts at 33 and 90 kpc are 
marked with vertical dotted and dashed lines, respectively. 



several objects which have been studied in detail ( Sato et al 



2008 2009). 



Trends in abundance ratios across cold fronts are an 
important missing piece of the puzzle in this ongoing debate. 
This is because sloshing effectively brings into contact gas 
which has been near the BCG and "pristine" gas at large 
radii, which has not been influenced (or has been influenced 
much less) by the BCG at late times. If the BCG has a 
large impact on the metal budget of the gas in its vicinity, 
this should be reflected in marked differences between the 
chemical enrichment patterns of the gas on either side of the 
cold front. To test this, we focus here on the northwestern 
cold front at 90 kpc, for which the good energy resolution of 
the Suzaku observations allows more precise determinations 
of the metal abundances than possible with other data sets 
(see previous work by e. g. Dupke & White 2003 Dupke 
et al.|2007| for the case of A496 with Chandra). 
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Table 1. Metal abundances on the bright and faint sides of the 
cold front. The Suzaku spectra extracted from the first 3 and last 
6 ellipses used to create the radial profiles were fit in parallel with 
the temperature and normalisations left free (3T/6T model) in 
order to determine the metal abundances on the bright and faint 
side of the cold front, respectively. The errors are given at the 
Ax 2 =l confidence level. 



Parameter 


bright side 3T 


faint side 6T 


O 


0.30 ±0.07 


0.17 ± 0.07 


Ne 


0.42 ±0.07 


0.31 ± 0.08 


Mg 


0.34 ± 0.07 


0.22 ±0.07 


Si 


0.36 ±0.04 


0.16 ±0.04 


S 


0.39 ±0.05 


0.23 ±0.05 


Ar 


0.55 ±0.10 


0.45 ±0.12 


Ca 


0.63 ±0.16 


0.60 ±0.18 


Fe 


0.41 ± 0.02 


0.25 ±0.02 


Ni 


0.35 ±0.14 


0.28 ±0.16 


X 2 /d.o.f. 


2730/2696 


3848/3678 



Since precise determinations of the metal abundances 
other than Fe require significantly more counts, we must 
substantially increase the sizes of the extraction regions. We 
thus can only determine one value of the O and Mg abun- 
dances for each side of the cold front, rather than being able 
to produce a radial profile. For this, we fitted in parallel the 
Suzaku spectra extracted from the first 3 and last 6 ellipses 
used to create the radial profiles in Fig. [3] to determine the 
metal abundances on the bright and faint side of the cold 
front, respectively. The temperature and normalisation of 
each ellipse spectrum was free in the fit, while all elemental 
abundances were coupled between the data sets in each case. 
This ensures that we account for variations in temperature 
across the large regions sampled by the spectra, which would 
otherwise bias the abundance measurements. Note that al- 
lowing the Fe abundance also to vary between the different 
elliptical extraction regions does not change the measured 
abundances for the other elements. The results are shown in 
Tabled 

We caution that the abundance of any element pro- 
duced predominantly by SNcc is unfortunately very diffi- 
cult to determine from the X-ray spectra. The Ne X line is 
blended with the strong Fe-L complex between 0.9-1.2 keV, 
which is very sensitive to the exact modeling of the tem- 
perature structure of the emitting gas. The O determina- 
tion is sensitive to the normalisation of the 0.2 keV galactic 
halo component, which is difficult to determine exactly be- 
cause the soft X-ray background around M87 is patchy. Rais- 
ing or lowering this normalisation by a factor of 2 changes 
the best-fitting O abundance on the faint side of the cold 
front by as much as 40%, while Mg changes only by 14% 
and the temperature by less than 0.02 keV. The Mg line 
is also blended with Fe-L emission and, although the Fe-L 
line emission around the Mg line is much weaker than that 
surrounding the Ne line, different spectral codes disagree 
on its modeling. The current collisional ionization equilib- 
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Figure 4. O/Fe and Mg/Fe abundance profiles. In black, values 
determined for several central annuli by |Matsushita et ah] ( |2003| ) 
with XMM-Newton. In blue, the results for the bright side and 
in red, for the faint side of the cold front. 



rium (CIE) plasma model implemented in SPEX (Kaastra 



et al.|1996 1 gives Mg abundances lower by a factor of 2 com- 
pared to the vapec model, which provides a better fit to the 
Mg/Fe-L line blend (e.g. |Matsushita et al.|[20"03"| , motivat- 
ing our choice of the latter plasma model for the present 
work. An improved version of SPEX containing an updated 
Fe-L line library, which is yet to be released, provides better 
agreement with vapec (J. de Plaa, private communication; 
see Table [2j|. The temperature structure moreover plays an 
important role in determining the elemental abundances, in- 
cluding Mg. If the metal abundances are obtained by fitting 
the integrated spectra of the regions in front of and beyond 
the cold front with single temperature models, there is an 
indication of a relative increase in Mg/Fe on the faint side 
compared to the bright side of the surface brightness jump, 
despite the systematic offsets in determining the absolute 
Mg/Fe between the different plasma models. However, ac- 
counting for the temperature structure by allowing the tem- 
perature to vary across several rings with fixed abundances 
shows that this trend is only an artefact. 

The results presented in Table [I] show a decrease both 
in the Fe abundance and in the abundances of all elements 
predominantly produced by SNcc across the cold front. Al- 
though the individual significances of the changes in the O, 
Ne and Mg abundances are small due to the large statistical 
uncertainties, the fact that all three abundances decrease 
by a similar relative amount (26-43%), which is approxi- 
mately the same as the relative change in the Fe abundance 
(40%), is a convincing indication that the measured differ- 
ences are reliable. Adding in quadrature the significances of 
the variations in O, Ne and Mg across the cold front sug- 
gests a change in the abundance of SNcc products at the 
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Table 2. Systematics in determining the Mg/Fe ratio. The 3T/6T 
model refers to fitting the first 3 and last 6 ellipses used to create 
the radial profiles in parallel with the temperature and normali- 
sations left free. 



Mg/Fe 


bright side 


faint side 


vapec 3T/6T 


0.81 ± 0.17 


0.88 ±0.29 


vapec IT 


0.77 ± 0.15 


1.14 ± 0.32 


spex IT 


0.42 ±0.15 


0.68 ± 0.32 


new spex IT 


0.97 ± 0.18 


1.28 ±0.40 



90% confidence level. This change requires the presence of a 
centrally peaked distribution not only for Fe but also for O, 
Ne, and Mg, and rules out models in which SNcc products 
are uniformly distributed throughout the ICM, while only 
SN la products show a radial gradient. 

Because all elemental abundances decrease by very sim- 
ilar relative amounts, we find no significant change in either 
Mg/Fe, Ne/Fe or O/Fe across the cold front. Comparing the 
Mg/Fe and O/Fe ratios on either side of the cold front with 
those measured at smaller radii by |Matsushita et aL| ( |2003[ ) 
using XMM-Newton data (also using the vapec model), we 
find a very slow increase across the entire radial range. The 
data are consistent with no increase in Mg/Fe from the very 
centre of the galaxy out to roughly 210 kpc (45'), as shown 



in Fig. Ul |Doherty et al. ( 2009 1 show that the stellar halo of 



M87 truncates at an average radius of ~150 kpc, therefore 
we are clearly probing the chemical enrichment history in re- 
gions both inside and beyond the sphere of influence of the 
galaxy. The very low Fe abundance outside the cold front 
(beyond 20' in Fig. [3| also strongly suggests that we are 
probing the chemical enrichment of the abundance floor at 
which the ICM metallicity is expected to flatten out at large 
radii. The roughly constant Mg/Fe and O/Fe abundance 
ratios indicate similar chemical compositions and, conse- 
quently, enrichment histories in all these regions, although 
within the large statistical errors we cannot rule out a recent 
influence of SN la in M87 on the metallicity budget of the 
Virgo ICM. 

Assuming the core-collapse model yields of Nomoto 



et al. ( 2006 I with solar metallicity of the progenitor, the 
best-fitting Mg/Fe ratio in the centre of M87 (based on the 
average of the three data points from Matsushita et al. 12003) 



corresponds to 27% of all supernova explosions contributing 
to the chemical enrichment being SN la, while the best- 
fitting Mg/Fe ratio on the faint side of the cold front suggests 
a relative contribution of TVsn i a /(iVsN la + Nsn cc ) ~20% 
(independent of the assumed SN la yield model). An enrich- 
ment of the gas on the faint side of the cold front (beyond 
90 kpc) purely by SNcc without any contributions by SN la 
would imply a Mg/Fe ratio of 2.6 solar and is ruled out at 
the 5.8<r level. 

The total mass of Fe contained within a 90 kpc ra- 
dius from the centre of M87 is approximately 2.3 x 10 s Mq 
( Bohringer et al.|2004| |. If 27% of all supernovae in the core 
region are SN la (based on the Mg/Fe ratio discussed above) 
and each SN la produces ~0.8 Mq of Fe (Iwamoto et al 



1999 1, while the rest are core collapse supernovae (which, for 
a solar metallicity of the progenitor and assuming a Salpeter 
IMF, produce appr oximately 0.09 Mq of Fe per explosion, 
Nomoto et al.|2006| ), then 2.2 x 10 s SN la and 6 x 10 8 SN C c 



are needed to produce all the metals observed in this region. 
For a Salpeter initial mass function (IMF) between 0.1 and 
100 M©, approximately one star with M > 8 M© which 
will explode as a SNcc is expected per 100 solar masses of 
stars formed. For a B-band mass to light ratio of around 



8 (Kronawitter et al 
10 k 



20001 and a blue-band luminosity of 



Lbq corresponding to the absolute B magnitude of M87 



(19901, the stellar mass 



of -22.14 cited by Peletier et al 

of M87 is 8 x 10 11 M Q and 8 x 10 y SNcc would have ex- 
ploded as the galaxy was forming, 13 times more than the 
6 x 10 8 SNcc needed to produce the metals in the central 
90 kpc from the M87 centre. The central Mg peak can thus 
be explained if only less than 10% of the SNcc products 
from the galaxy formation phase are retained rather than 
being uniformly distributed into a flat abundance profile. 
The expected present SN la rates in elliptical galaxies are 
still relative ly uncertain, from 0.18 ± 0. 06 SNU ((C appellaro 
19991) to 0.3 6j~^l ± 0.02 SNU ([Sh aron et al.f|2007) 



et al. 



or 0.28T^ SNU (jMannucci et al. 2008}, where 1 SNU 



equal to one supernova explosion per century per 10 10 Lbq 
of the galaxy. Since redshift z=l, corresponding to a look- 
back time of 7.7 Gyr, between 10 s and 4.6 x 10 s SN la could 
have enriched the centre of M87 within the current uncer- 
tainties, between half to twice the amount needed in total. 
It is therefore possible that a fraction of the SN la prod- 
ucts in the central peak stem from the proto-cluster phase 
rather than recent SN la in M87, although that requires the 
SN la rates to be on the lower end of the measured confi- 
dence intervals and not to have been larger in the past, as 



predicted by e.g. Renzini et al. (19931; using 0.28 SNU, the 



(2008), exactly 



most recent measurement of Mannucci et al 
the required number of SN la, 2.2 x 10 s , would explode since 
z=l, thus all the Fe would be produced only by the central 
galaxy at late times. Another possibility to explain the con- 
stant Mg/Fe ratio in this case is that some of the Fe outside 
the metal peak was produced recently by SN la explosions 
of intra-cluster stars. For a diffuse intra-cluster lig ht with 
an average apparent magnitude of 27 mag/arcsec 2 (Mihos 
et al.||2005| |Kormendy et al.||2009|, 30% of the Fe outside 



the 90 kpc cold front and within the radius covered by the 
XMM-Newton mosaic could have been produced by SN la 
since z=l, for a SN la rate of 0.28 SNU. 

We plot in Fig.lBlthc measured Fe/S vs. Fe/Mg ratios for 
the M87 centre ( Matsushita et aT]|2003 1 and on either side 
of the observed cold front using our measured values. The 
Fe/Mg, as discussed above, is an important diagnostic of the 
relative contribution of SN la to SNcc • Using this ratio and 
assuming the same core-collapse model yields as above, we 
can determine the amount of Si-group elements produced by 
the SNcc and consequently the ratio of Si-group elements 
to Fe produced by SN la, which is an important diagnostic 
for the explosion mechanism of the latter. We overplot in 
Fig. [5] three different abundance ratio models based on the 
deflagration ( W7) and delayed deton ation (WDD1, WDD2) 
SN la yields ofllwamoto et al. ( 1999 1. We used S as a repre- 
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Figure 5. Fe/S vs. Fe/Mg abundance ratios on the two different 
sides of the cold front (this work) and in the centre of M87 | |Mat-| 
|sushita et al,|2003| | . Overplotted are supernova enrichment models 
obtained by combining various relative numbers of SN la (using 
the WDD1, WDD2 and W7 yield models of |Iwamoto et al.|1999| 
and SNrjc (whose yields are modeled based on |Nomoto et ah] 
|2006| assuming a solar metallicity of the progenitor and averaged 
over the Salpeter IMF). 



sentative element for the Si-group because the Si abundance 
is more easily biased by calibration issues around the Si edge 
and by modeling the underlying Si instrumental line. We 
find marginal evidence supporting the scenario proposed by 
12002), in that the data points at larger 



Finoguenov et al. 



radii favor a W7 model whereas the abundance ratios at the 
centre of M87 favor the WDD1/WDD2 SN la yields. 



5 CONCLUSIONS 

Two cold fronts are observed in M87: one towards the SE 
at a radius of ~33 kpc, and one towards the NW with a 
radius of ~ 90 kpc. The opposite and staggered placement 
of these cold fronts, as well as the absence of a remnant 
subcluster, make sloshing the most viable explanation for 
the presence of these cold fronts. The ICM distribution in 
and around M87 is very similar to results obtained from nu- 



merical simulations of sloshing in cool core clusters (Tittley 



& Henriksen 2005 Ascasibar & Markevitch 2006). For the 



mass profile of M87, the time needed for two fronts at these 
radii to become out of phase by 180 degrees is approximately 
0.4 Gyr. 

We present accurate projected temperature profiles 
which show that the temperature is higher on the X-ray 
fainter side of each cold front, opposite from what would 
be expected in the case of shock fronts. We also show that 
the metallicity drops sharply outside each cold front. Such a 



discontinuity in the metal abundance profile is expected be- 
cause sloshing displaces the central, cooler, more metal-rich 
gas and brings it into contact with hotter, more metal-poor 
ICM in the outskirts. In this way, sloshing also contributes to 
transporting the heavy elements from the centre outwards. 

We estimate that the inner cold front transports 4.5 x 
10 6 Mq of Fe, while for the outer cold front the total mass 
of Fe transported amounts to 13.7 x 10 6 Mq. These masses 
represent in each case the same fraction (~6-8%) of the total 
Fe mass available within the radial range affected by slosh- 
ing, which is much larger than the amount of Fe transported 
by AGN-related processes, amounting to only 1.5 x 10 6 Mq 
of Fe. 

We find no strong indications for a change in the Mg/Fe 
and O/Fe ratios across a large radial range including the 
outer cold front, given a proper modeling of the tempera- 
ture distribution within the extraction regions. This requires 
the presence of a centrally peaked distribution not only for 
Fe but also for O, Ne, and Mg, and rules out models in 
which SNcc products are uniformly distributed throughout 
the ICM, while only SN la products show a radial gradient. 
Since stellar mass loss in the central galaxy typically con- 
tributes only small amounts of heavy elements to the ICM 
metal budget, the peaks in the distribution of SNcc prod- 
ucts are likely to stem from star formation, triggered as the 
BCG assembled during the protocluster phase, as proposed 
by Simionescu et al. (20091. The central Mg peak can be 



explained if only less than 10% of the SNcc products from 
the galaxy formation phase are retained rather than being 
uniformly distributed into a flat abundance profile. 

The very low Fe abundance of only 0.25 ± 0.02 solar 
(0.17 ± 0.01 relative to |Anders fe Grevessep989[ ) suggests 
that the gas outside the cold front is representative of the 
pristine outskirts of typical clusters. This means that we 
are currently observing the very mechanism which spreads 
the metals into these pristine regions, broadening the metal 
distribution to that usually seen in relaxed clusters. More- 
over, we are able to determine the chemical composition of 
this pristine gas, and find that probably around 20% of all 
supernovae contributing to its enrichment are SN la, that 
the metal abundance ratios of this gas are very similar to 
those near the centre of the BCG, and that an enrichment 
by SNcc products only can be excluded with a very high 
significance. 
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